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ABSTRACT 


Development of flats on the wheels of rail*? md wa « ons is 0ne ° f ^ maj ° r 

„ . , t • r , tion . Wheel flats are flat zones on the wheel 

causes of track and ride quality detenorau^ 

, , . . . , .... f t u e wheel on the rail when the brakes lock. A 

tread caused by unintentional sliding or rn c 

fresh wheel-flat has the shape of a geometnc chord a. Are wheel circumference. Forces 

from the wheels with flats can be up to four i*» *** ,han those ftom n0rma ' ^ 

. ,1 wheel-flat and rail interaction dynamics 

The objective of the present work is to mou^ 

, , _ . , ... , tn understand experimentally observed response 

through a finite element model, in order to u 

. . procedures are validated for free and forced 

of an instrumented railway track, lne p iu 

. . , • ntl with multi-span continuous beam analytical 

vibration characteristics through comparison 

ar e generated and response is computed for 

models. Typical wheel flat loading patterns 

, . . . . fl ... cjnports and damping. Response patterns are also 

various cases of rigid and flexible rail supp^ 

compared with test results from an instrument track. 
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CHAPTER 1 


INTRODUCTION 

Heavier trains and higher speeds accentuate the detrimental influence of imperfections in 
vehicles and ti'acks on the dynamic interaction between the two. Large dynamic load may 
be excited causing significant damage to both track and vehicle. Common imperfections 
are irregularities on the running surface of rail and wheel and deficiency in the support of 
track structure. Isolated defects such as wheel flats or badly aligned rail joints result in 
high frequency vibration, while for example loss of ballast under a single sleeper mostly 
affects the low frequency response of the track. The dominating frequency excited by a 
periodic irregularity, such as a sinusoidal corrugation of the railhead or the wheel tread, 
depends on the wavelength of the irregularity and the speed of train. 

The influence of out of roundness of a wheel on the dynamic train/track is often 
considerable. A severe type of out of roundness is the wheel flat, which is a flat zone on 
the wheel tread caused by unintentional sliding of the wheel on the rail when the brake 
locks. A fresh wheel-flat has the shape of a geometric chord at the wheel circumference. 
However, when the brake unlocks the comers of the flat profile soon become rounded. 
The forces from the wheel flats are considered to be the decisive contributory factor 
affecting track degradation and life. Wheel forces primarily influence the rails and 
concrete sleepers. However, it is known that forces pass through the entire substructure 
down to the subsoil. 

The objective of the present work is to model, the wheel-flat and rail interaction 
dynamics through a finite element model, in order to understand experimentally observed 
response of an instrumented railway track. The finite element model is developed using 
commercially available software (NISA). The finite element procedures are validated for 
free and forced vibration characteristics through comparison with multi-span continuous 
beam analytical models. Typical wheel flat loading patterns are generated and response is 
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computed for various cases of rigid and flexible rail supports and damping. Response 
patterns are also compared with test results from an instrumented track. 

The general problem on rail-wheel interaction attracted attention as early as 1 847, when a 
British commission was appointed to study dynamic effects in railway bridges. The 
differential equation of deflection was established by neglecting inertia of the bridge and 
by considering the moving load as a concentrated mass. An exact solution of the equation 
was obtained by Stokes (1849). The next step was made by Timeshenko (1922), who 
solved the problem of a pulsating forces moving with uniform velocity along a beam. 
Other work has been done, including that of Lowan (1935), who studied the case in 
which the velocity of traversing force is not constant. The foregoing investigations were 
all limited to single span structures. After that Ayre, Ford, and Jacobsen (1950) have 
studied transverse vibration of a two - span beam under action of a moving constant 
force. Some of the more recent works include work done by Cai, Cheung and Chan 
(1988) who studied dynamic response of infinite continuous beam subjected to moving 
force. They made the assumption that slope and moment at right end of last span is equal 
to slope and moment at left end of first span respectively. For multi-span beams most 
researchers have used Finite Element methods. Nielsen and Igeland (1995) has 
investigated vertical dynamic behaviour for a railway bogie moving on a rail by use of an 
extended state-space vector approach in conjunction with a complex modal superposition 
for the track. Effects of imperfections on the running surface of wheel and rail were 
studied by assigning irregularity function to these surface and calculated contact forces 
for moderate vehicle speed in the presence of wheel flat. Thambiratnam and Zhuge 
(1996) have developed simple finite element model to analyze simply supported beams 
on an elastic foundation of any length, where the foundation has been modeled by springs 
of variable stiffness. The effect of the speed of moving load, foundation stiffness and 
length of beam on the response of the beam have been studied and dynamic amplification 
of deflection and stress have been evaluated. Abu Hilal and Zibdeh (2000) have 
investigated vibration of elastic homogeneous isotropic beam with general boundary 
condition traversed by moving force. They however, considered only single span beams. 
Wu and Shih (2000) have studied dynamic response of railway and carriage under the 
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high speed moving loads. Dynamic response of the rail and carriage due to action of 
multi-roller carriage was determined by means of the finite element method. Wu, and 
Thompson (2001) have investigated vibration analysis of railway track with multiple 
wheel on the rail. However the focus is on railway rolling noise. 

A wheel flat is caused by the brake or some other mechanism, interfering with the normal 
angular velocity of the wheel. When brakes do not function properly, the wheel ceases to 
rotate and a flat develops due to plastic deformation. The dimensions of the wheel flat, its 
form and position on the tread depends on length of the wheel skid, axle loading, running 
speed, shape of rail profile, position of the wheel on rail, wheel material and the material 
coefficient of adhesion between the wheel and the rail material. 

Wheel flats cause large dynamic forces to the track. According to popular practices, flats 
on wheels should be restricted to a length of 60 mm and a depth of 0.9-1. 4 mm. In this 
range dynamic wheel load increases by approximately 30 kN/mm for timber sleepers and 
50 kN/mm for concrete sleepers. Wheel flats are also known to cause very high 
acceleration (approx. 500 g) on the track. These high acceleration become input to the 
rolling stock, which will induce greater forces than permitted and damage the suspension 
system frame and the body of the rolling stock. These forces may also cause bearing 
damage and hot axles boxes. Wheel with flats causes local unevenness on the tread. The 
foreign bodies may penetrate in to the tread due to wheel deformation. Deformed wheels 
are also damaged by the heavy impacts. In addition, excessive vibration and unpleasant 
air borne noise is caused which is irritating to the passengers. 

In the present work the rail track has been modeled as a multi-span continuous beam with 
rigid or flexible supports at regular intervals. The free vibration analytical model is 
described in Chapter 2 along with the finite element model. Forced vibration under a 
single moving load is discussed in Chapter 3 and the finite element procedure is validated 
by comparison with results from an analytical solution. Rail response under a moving 
train, with representative data, is obtained from the finite element model in Chapter 4, for 
cases of wheel with flat and without flat. The results are compared with those from a test 
site. Conclusions and scope for future work are given in Chapter 5. 
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CHAPTER 2 


FREE VIBRATIONS OF 
MULTI-SPAN CONTINUOUS BEAMS 

Mathematical models of varying degrees of complexity have been developed by 
researchers to analyze rail-wheel interaction. The focus of the present study is to 
develop a Finite Element model of the rail in order to study its response under typical 
wheel-flat loading. The free-vibration characteristics of rails are functions of its 
geometric and material characteristics. It is customary to model the rail as a 
continuous beam with multiple supports. These supports to the rail come in the form of 
wooden or concrete sleepers, which are laid out periodically across the rail (Fig.2.1). 
The sleepers can be modeled as rigid supports or flexible springs. Free vibration 
characteristics of rails modeled as multi-span continuous beams are discussed in this 
chapter. The analytical formulation is based on the procedure suggested by 
Timoshenko (1955) and is employed to validate the results from the finite element 
model. 

2.1 Governing Equation 

Consider the case of a continuous Euler - Bernoulli beam with N spans simply 

supported at the ends and (jV- 1) intermediate support. Let lij 2 In be the lengths of 

consecutive spans, the flexural rigidity of the beam being the same for all spans. 
Taking the origin of coordinates at the left end of each span as shown in Fig. 2.2, 

the equation for transverse vibration for each span is given by 

El 0, 0 <*,</,, r = 1,2., N (2.1) 

dx r dt 

Suffix r denotes the rth span, El, p and A denote, respectively, the flexural rigidity, 
mass density and the cross-sectional area, w is the transverse deflection of each span, 
x is the local coordinate along the axis of the rth span and t is the time . The continuity 
and equilibrium condition at the intermediate support points of the continuous beam 
require following relations: 



w r (x r =0,0 = 0, 
w r (x r = /,.,?) = 0, 


, 5w„,. 

— (x r =l r ,t) — — (x r+1 =0,/), 

C/X r <7X r +l 

d 2 w r d 2 w +1 

~T ~ ( X r = hd) = 2 ( X r+1 = 0,0, ^ = 1,2, JV-1 (2.2) 

3x r 5x0+1 

2.2 Response 

Let the solution of equation (2.1) be given by 

w r = X r (x) * (5, sin cot + B 2 coscot ) (2.3) 

where 2C- is normal function for span r and co is the natural frequency and B ; and B? 
are constants. 


The normal function for the span r will be 

X r = a r coskx + b r coshAx + c r sin Ax + d r sinhAx (2.4) 

Substitution of the first of boundary conditions given in (2.2) gives 
a r = - b r , 
so that 

X r = a r (cosAx - coshAx) + c r sinAx+<7 r sinhAx (2.5) 

In the above a r , c r , d r remain constants to be determined from the condition at the 
ends of the span and 



The consecutive derivatives of (2.4) will be 

f 

X r = a r A{sinAx + sinhAx) + cy A cos Ax + coshAx (2.7) 

X" =-a r k 2 (coskx + coshkx)-c r k 2 sinkx+d r k 2 sinhAx (2.8) 

From the condition at the simply supported ends of the beam it can now be concluded 
that 





ai-aN+i-O. 


Considering the condition at the right end of span r we have, 

),=/,. = 0 > 

(x;.) x=lr =(X' r+l ) x=0 , 

{X" r ) x=lr =(^; + 1 ), = o. 

Now from equations. (2.4), (2.7), (2.8) one gets 
a,, (cos kl r - cosh kl r ) + c r sin kl r + d r sinh kl r = 0, 

- a r (sin kl r + sinh kl r ) + c r cos kl r + d r cosh kl r - c,. +l + d r+i , 
a r (cos kl r + cosh kl r ) + c r sin kl r - d r sinh kl r = 2 a r+] , 


Equations (2.10) - (2.12) can be solved to obtain 

_ a r+l - cos U r 
sin kl r 


■ a r+] + a r cosh kl r 
sinh kl r 


(for sin kl r not zero) 


From the above 

c r + d r = a r (coth kl r - cot kl r ) - a r+] (cos ech kl r - cos eckl r ) 

Using notations: 

coth kl r - cot kl r = (p r 

cosec hkl r - cosec kl r = y/ r 

we obtain 

c r +d r = a r (j> r -a r+x y r 


Similarly, for the span r+1, 
c r+l +d r+ 1 =a r+ \<f ) r+ 1 -a r+2 f r +\ 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 

(2.15) 

....(2.16) 


Substituting equations (2.13) and (2.16) in equations. (2.1 1), we obtain 



Q-yty r ®r+\($r ^r+1 ) ^r+lWr+X — ^ 


(2.17) 


Writing analogous equations for each intermediate support, the following system of 
(N-J) equations will be obtained: 

~a 2 (0 1 + <f> 2 ) + a 3 i// 2 =0. 

a 2 y/ 2 - a 3 (<p 2 +<Pz) + = 0, 


~ q n( ( Pn-\ + </>\ ) - 0, (2.18) 

The frequency equation for the free vibration will be obtained by putting the 
determinant of the above equations equal to zero. 

If the length of each span is same means UJ. 2 =lu = l then the relation for first natural 
frequency is 
k\l - re 

Natural frequency and mode shapes can be found by solving equations (2.18), (2.6) 
and (2.5). 

2.3 Finite Element Model 

The analytical formulation described in the previous section is employed to compare 

the results from a finite element model for a multi-span beam. The finite element 

algorithm is to be used further to model the rail track. A simple 2-D beam-element, 

with eight elements per span, has been employed in the finite element model. There 

are two nodes per element and the element has 3 degrees of freedom namely, 

u x - displacement in X-direction, 

u y = displacement in Y-direction, 

r z = rotation about Z-direction. 

u x = 0 at all nodes and u y =0 at supports 



Two beams with 15 and 25 spans respectively, as shown in in Figs.2.3.(a) and (b) 
were chosen to carry out the free vibration modeling exercise and comparison of 
results with the analytical model described in sections 2.1 - 2.2. For Beam 1, with 15 
spans the total number of elements was 120 (121 nodes), while for Beam 2, with 25 
spans the total number of elements was 200 (201 nodes). 

2.4 Numerical Results 

The numerical exercise was carried out for the following set of data 
Young modulus of beam E= 2 x 10 n N/m 
Cross section Area of beam A = 9 x 1CT 4 m 2 

Area moment of inertia I = 6.75x10 m 
Density of beam material p = 8000 Kg/rn 
Length of each span / = 1 m 


The eigenvalue problem of equations (2.18) is solved to obtain natural frequencies 
and mode shapes of beams with 15 and 25 spans. Eigenvalue analysis of the beams 
was carried out in NISA using conventional subspace iteration and consistent mass 
formulation. First ten natural frequencies thus obtained analytically are given for 15 
span and 25 span beams in Table 2.1 and Table 2.2 respectively. The results obtained 
from the FE model developed using the commercially available software NISA are 
also given in the tables. The close agreement between the two sets of results validate 
the finite element modeling procedure in NISA. The corresponding mode shapes for 
the two beams are shown in Figs. 2.4 (a), (b) and 2.5 (a), (b) respectively. 



Table 2.1 Natural frequencies for 15 span beam 


Mode No. 

Analytical model 
(rad/s) 

FE model 
(rad/s) 

1 

427.37 

427.22 

2 

432.81 

432.65 

-> 

J 

448.74 

448.58 

4 

474.14 

473.97 

5 

507.62 


6 

547.68 

547.47 

7 

592.91 

592.68 

8 

642.01 

641.76 

9 

693.74 

693.47 

10 

746.81 

746.51 


Table 2.2 Natural frequencies for 25 span beam 


Mode No 

Analytical 
model (rad/s) 

FE model 
(rad/s) 

1 

427.37 

427.21 

2 

429.33 

429.18 

J 

435.18 

435.02 

4 

444.76 

444.60 

5 

457.84 

457.67 

6 

474.14 

473.97 

7 

493.35 

493.17 

8 

515.15 

514.96 

9 

539.21 

539.00 

10 

565.22 

565.001 







Fig 2.1 Rails supported on sleeper pads 
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Fig 2.2 Multi-span continuous beam 
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Fig 23 (a) Finite element model of 15 span beam 
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Fig 2.3 (b) Finite element model of 25 span beam 
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Vertical modal displacement 




Fig 2.4 (a) Mode Shapes for 1 5 span beam - Analytical formulation 
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CHAPTER 3 


VIBRATION OF MULTI-SPAN CONTINUOUS BEAMS 
UNDER TRAVELLING LOAD 

The response ot multi-span continuous beams subjected to moving loads is discussed in 
this Chapter. Using principle ot virtual work general expressions for forced vibration 
response (1 imoshenko, 1 955) are presented. These are consequently employed to express 
beam vibrations caused due to constant magnitude moving loads. The analytically 
obtained response is compared with that from the FE model. 

3.1 Governing Equation 

Representing the lateral vibration of a beam ( Fig 3.1 ) subjected to external disturbing 

forces, by a series 

vp = /q A i + jUn A 2 + P^X ^ + (3-1) 

in which A'/, AT... are normal functions of the beam correspondent to the end condition of 
the entire vibrating beam and Py,Mi. are time functions. 


P 



/ 


1 I — >1 1 i 1 1 1 1 1 

Fig 3. 1 Multi spans continuos beam under a moving force 


The inertia forces per unit length of the beam can be expressed as 

- pAw = -pA 'Zp m X m ( 3 - 2 ) 

For a virtual displacement 

x m (3 - 3) 
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the virtual work done by inertia forces is 


P A l w Swdx =-pA }u m SjU m \x 2 m dx 


(3-4) 


where L is the total length of beam. 


The strain energy of bending beam at any instance is 

FI 1 tl 2 w 

L )^' 

- o dx~ 

which can be expressed as 
FI 1 L T 

=4- Z XI- dx 

- m () 


Similarly, the virtual work done by elasticity forces can be expressed as 


dK 


dju 


dp 


m 


FI L -J 


(3.5) 


m 


For the disturbing force P applied at the distance c from the left end of beam of first span, 

the virtual work is 

(3-6) 


Applying principle of virtual work to expressions (3.4), (3.5), (3.6), the equation of 

motion is obtained as, 

I. 

X".~dx 

P 


Mm "*■ 


r/ \Krdx 
EI o 


pA l : v 


JU 


m i 




(3.7) 


\Xjdx pA\X m dx 

0 0 


Employing the notation 
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(3.8) 



EJ_ 

pA 


1 9 

\X"Jdx 
o 

\x,;-dx 


0 


equation (3.7) can be written as 
..2 P 

Mm +A m ~Mm = - y— — (-AT m ) x ~ c 

pAjX m ~dx 

o 

The genera! solution of equation (3.9) is 

p t 

Mm ~~ ^ m ^in ~ J J (X m ) sin X m (t — q )<7/| 

0 


(3.9) 


(3.10) 


where and D m are constants, which can be found from initial conditions. The third 
term represents the vibration produced by disturbing force. 


For zero initial conditions, we get 

M m — - J P( X m ) sin X m (t - 1, )dt x (3.11) 

pAX m J X m 2 dx 0 

o 


3.2 Response to Moving Force 

Consider the disturbing force P to be moving along the beam at velocity V. Let the force 
be at the left support of the first span of the beam, at the initial moment (t = 0). Then at 
any other time instant moment r = t u the distance of the force from the left support of the 

first span will be 
c '----Vt, 

The virtual work done by the moving force 

P^(X m )„ n , 0 . 12 ) 

Through comparison of equations (3.12) and (3.6) and substituting (X m ) Y=K/| for 
{X n ) xmc in equation (3.1 1 ), the beam response due to the moving force can be written as 
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Mm 


pAX 


,nU„ 

0 


l dx 


1 m ) sill X m (t — /■] )<S?f, 


x~Vt\ 


(3.13) 


3.3 Numerical Examples 

The response is numerically computed through the procedure outlined and compared with 
the hinite h lenient Model developed in NISA. This comparison is carried out to validate 
the finite element model, hoi the multi-span beam, earlier specified in section 2.4, the 
forced vibration response is numerically computed for the following set of data: 


Case 1: 

Total number of spans N 15 

Moving load P ION 

Total number of modes considered 12 

Velocity. I ' of the moving load 2.5 m/s 


The response at the mid point of the central span (r = 8) obtained through the analytical 
procedure and the finite element method (NISA) are plotted in Figs 3.2 (a) and (b), 

respectively. 


Case,2: 


Total number of spans N 15 

Moving load P 10 N 

Total number of modes considered 12 

Velocity, I ' of the moving load 5.0 m/s 


The response at the mid point of the central span (r = 8) obtained through the analytical 
procedure and the finite element method (NISA) are plotted in Figs 3.3 (a) and (b) 5 

respectively. 
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The response obtained from the analytical formulation and the FE models can be seen to 
be identical toi both cases. This is taken as a pointer to the correctness of the finite 
element modeling procedures developed using NISA. 
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Fig 3.3 (a) Response for Case 2 - Analytical method 
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CHAPTER 4 


RAIL RESPONSE UNDER WHEEL-FLAT LOADING 

Rail response under wheel-flat excitation has been investigated. Initially the rail vibration 
under loading due to normal wheels is investigated. The analysis is consequently 
extended to incorporate loading due to wheels, which have developed flats. Modeling is 
done through the Finite Element package NISA and the wheel loads are treated as 
moving loads. 


4.1 Railway Data 

Analysis is earned out for typical data obtained from Indian Railways. A railway train 
comprising of an engine and one wagon are considered. A schematic of the train is given 
in Fig.4.1 (a). The engine has 12 wheels (six wheels on either side). On each side there 
are three front wheels and three rear wheels. Each wagon has eight wheels, four on either 
side, each side having two front wheels and two rear wheels. Analysis is carried out for 
only one side of the railway track (the response on the other side being identical). The rail 
cross-section is shown in Fig. 4.1 (b). The data employed for the analysis is listed below: 


Cross section Area of rail 

66.15 x 1CT 4 m 2 

Moment of inertia about N.A: 

2158x 10' 8 m 4 

Distance between two-sleeper pad 

0.6 m 

Young modulus of rail 

2.0678 x 10 n N/m : 

Density of rail 

7600 Kg/m 3 

Speed of train 

36 km/h (=T0m/s) 

Wheel diameter 

1 m 

Total load on each wheel of engine 

92214 N 

Total load on each wheel of wagon (empty) 

29062.125 N 

Total Number of wheels in the engine 

12 

Total Number of wheels in the wagon 

8 

Equivalent damping ratio of rail (rigid support) 

0.05 

Stifness of sleeper pad 

2.5 x 10 8 N/m 

Equivalent damping ratio of rail(flexible support) 

0.1 
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The weight of the engine and wagon is taken to act on the rail through the wheels as point 
forces. The respective weights of the engine and wagon are assumed to be equally 
distributed among the wheels. At the initial time instant the first wheel of the engine is 
assumed to be at node 1 of span number 1 of the rail. Table 4.1 gives the wheel loading 
data for all the ten wheels ( 6 from engine and 4 from the wagon). 


Table 4.1 Wheel loads, distances and arrival times at node 1 

(Train velocity =10 m/s) 


Wheel 

No. 

Arrival time 
at node 1 (s) 

Wheel load 
(*N) 

Distance from first 
engine wheel (m) 

1 

0.00 

92.214 

0.0 

2 

0.21 

92.214 

2.1 

o 

J 

0.42 

92.214 

4.2 

4 

1.05 

92.214 

10.5 

5 

1.26 

92.214 

12.6 

6 

1.47 

92.214 

14.7 

7 

1.79 

29.062 

17.9 

8 

1.99 

29.062 

19.9 

9 


29.062 

24.4 

10 


29.062 

26.4 


Total time taken to cross rail by train = (rail length + train length) /velocity 

= (9 + 26.4)710 = 3.54 s 


4.2 Finite element modeling 

Two models of rails are developed, namely (i) rail with rigid supports (ii) rail with 
flexible supports. A rail track of 15 spans is considered. The length of each span is 0.6 m. 
Different scenarios with wheel flat falling on the supports and in-between supports, are 
discussed. Details of the FE model are given below: 


(i) Rail with rigid support 

Type of element 2-D beam 

Order of element: 1 (linear) 

Number of degrees of freedom per node of beam element 3 (u x , u y , r 
Number of nodes per element: 2 

Total number of elements in rail 1 1 8 
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Total number of nodes in rail 
Total length of rail 

u x = 0 is taken at all nodes of rail and u y = 0 at supports only 


119 
9 m 


(ii) Rail with flexible supports 
(in addition to the data of case (i) ) 

Type of element of spring 2-D translation 

Order of spring element 1 (linear ) 

Number of degree of freedom per node of spring element: 2 (u x , u y ) 

Number of nodes per spring element 2 

Total number of spring elements 15 

u x =0 is taken at all nodes of rail and u y =0 at one end of spring element only 

Finite element modeling of rail with rigid supports and flexible supports are shown in Fig 

4.2 (a) and Fig 4.2 (b) respectively. 

4.3 Free vibration characteristics of rail 

Eigenvalue analysis of the rail was carried out in NISA using conventional subspace 
iteration and consistent mass formulation. Table 4.2 gives the first five natural 
frequencies of rail for the rigid support and flexible support models. The corresponding 
mode shapes are given in Figs. 4.3 (a) and (b). 

Table 4.2 First five natural frequencies of rail (rad/s) 


Mode No. 

Rigid support 

Flexible support 

1 

7842.73 

2854.00 

2 

7941.44 

2854.11 

o 

J 

8230.34 

2862.73 

4 

8690.14 

2891.54 

5 

9295.03 

2955.22 


4.4 Wheel flat loading 

The loading pattern caused due to out-of-roundness of the wheel depends on the severity 
and shape of the damage on the wheel. The out-of-roundness is however idealised by a 
flat on the surface, which would be seen as a chord on circular profile of the wheel. Such 
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a flat results in loss of contact between the wheel and the rail for a short time and then 
regaining of contact. This idealised lift and hit phenomenon is shown in Fig. 4.4. When 
there is no contact, the wheel rises above the rail until contact is recovered resulting large 
impact force. However, due to the bogie suspension system, the loss of contact between 
the wheel and the rail is not total as the wheel moves downwards and the rail upwards to 
compensate for the absence of wheel material. The overall effect is that the contact 
surface initially decreases and later due to inertia, the wheel continues downward motion 
resulting in an increase of contact force. The maximum and minimum values of 
normalized contact force between the wheel and the rail, for various speeds were 
calculated by Neilsen and Igeland (1995) and are reproduced here in Fig. 4.5. 

4.5 Vibration of rail under moving train 

Rail vibration under the moving train has been computed in terms of deflection and 
acceleration. The ten point loads from the ten wheels are moving on rail simultaneously. 
Response has been computed for the following cases: 

1 Undamped System 

(a) without flat 

(b) with flat 

2 Damped System 

(a) without flat 

(b) with flat 

The wheel flat is assumed to be present on the 1 0 th wheel of the train. The chord length of 
the flat is taken as 60 mm. The forced vibration analysis is carried out in NISA and 30 
modes are taken under consideration. Modes beyond 30 are found to have negligible 
influence on the response for the train velocity of lOm/s (=36 Km/hr), under 
consideration. 

The wheel-flat hits on the rail on the 5 th and then on the 10 th spans. In order to improve 
accuracy of results 20 element are taken on these spans, while on the remaining spans the 
number of finite beam elements are 6. Thus, the response is computed for 119 nodes on 
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the rail. The loading patterns for the case with and without wheel flats are shown in 
Figs.4.6 (a) and (b) for 5 th and 10 th span respectively. The flat wheel looses at node 33 
and regains it at node 35. Again, it looses contact at node 82 to regain it at node 84.. The 
nodal distances and arrival time of the 10 th wheel at these nodes along with the tenth node 
(which sees normal wheel behaviour) are listed in Table 4.3. 


Table 4.3 Nodal distances and arrival time of 10 th wheel 


Node 

Number 

Distance from 
Node 1 ( m) 

Arrival time of the 

10 th wheel (s) 

10 

0.90 

2.730 

33 

2.64 

2.904 

34 

2.67 

2.907 

35 

2.70 

2.910 

82 

5.79 

3.219 

83 

5.82 

3.222 

84 

5.85 

3.225 


Node 35, where the hit occurs for the first time, lies in the centre of a beam span (number 
5). Node 84, where the hit takes place for the second time is close to a support (i.e. a 
sleeper). The vibration of rail is expected to be more due to the first hit, in comparison to 
that due to the second hit. 


For the undamped system with rigid supports the displacement response at nodes 10, 33 
and 35 is shown in Figs. 4.7 (a), (b) and. (c) corresponds to the normal case (without flat) 
and Figs 4.8 (a), (b) and (c) gives the response in the case of a wheel with flat. 

The corresponding acceleration response is shown in Figs. 4.9 (a), (b) and (c) and 4.10 
(a), (b) and (c). It can be seen that while each wheel causes distinct deflection as it passes 
through each node, the acceleration response does not provide a clear picture since it 
comprises of a large number of frequency components. It can be said that measurement 
of displacement response through strain gauging techniques can provide better 
information than that obtained through measuring acceleration, using accelerometer 
sensors. 
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From the plots of Figs. 4.8 (b), the impact caused by the wheel flat is clearly visible. The 
maximum deflection due to the engine wheels is 47.82 /urn at node 35. Maximum 
deflection, at node 35, due to the wagon wheels is 15.59 //m in the absence of any flat, 
while it rises to 27.o5//m, when a flat is present. Occurrence of this increase in the 
maximum deflection can serve as an effective tool for wheel-flat detection. The zoomed 
picture of the response under wheel flat loading is depicted in Fig.4.11, which clearly 
indicates the lift and hit mechanism of the wheel flat discussed earlier. 

The response for the damped system is shown in Figs. 4.12 (a), (b) and (c). The inclusion 
of damping in the analysis, as can be seen, removes the free vibration components, from 
the response plots. The zoomed picture of the response under wheel flat loading is given 
in Fig.4.13. 


Computation is also carried out for cases where the rail support has been taken as 
flexible. The response in this case is shown in Figs. 4.14 (a), (b) and (c). As expected 
deflection is higher in this case in comparison to that of the rail on rigid supports. 

The response at nodes 82, 83, 84 (for rigid supports and with damping) is shown in Figs. 
4.15 (a), (b) and (c). Comparison with the response at nodes 33, 34, 35, discussed earlier 
reveals that the deflection levels are lower and the impact due to wheel flat is also less 
pronounced. This is due to the fact that while nodes 33,34,35 lie in the centre of a beam 
span, nodes 82, 83, 84 are closer to a support (i.e. sleeper). 

The maximum deflections caused by the wagon wheel, in various cases discussed above 
are listed in Tables 4.4 for quick comparison. 


Table 4.4 Maximum deflection due to wagon wheel 

(at node 35 ) in ju m 


Type of support 

Without flat 

With- 

flat 

Undamped 

Damped 

Undamped 

Damped 

Rigid support 

15.59 

14.89 

27.35 

26.26 

Flexible support 

221.26 

79.67 

203.00 

109.75 
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4.6 Remarks 

The finite element algorithm has been developed to understand the experimental results 
obtained from a test site. The tests have been carried out through strain-gauge 
instrumentation of railway track. Signals have been obtained from these sensors during 
passage of the train under discussion. Strain-gauge response obtained at the test site is 
shown in Fig. 4.16 (a) - for all normal wheels, (b) - for a 60 mm flat on wheel number 10 
Fig 4.16 (c) - magnified plot of the response in (b). The deliberately introduced flat on the 
test wheel has been shown in Fig 4.17. The test results help to show that the Finite 
Element algorithm has been correctly able to model the rail vibration phenomenon under 
wheel-flat excitation. 
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(All dimensions are in mm) 

Fig 4.1 (b) Rail Section 
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2nd Mode 


3rd Mode 


4th Mode 
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Fig 4.3 (b) Mode shapes of Rail -Flexible supports 
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(c) 

Fig 4.4 Lift and hit phenomenon of Flat wheel 

(a) Start of lift, (b) lift condition and (c) hit condition 
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0 20 40 60 80 100 120 


Speed (Km/h) 

(b) 

Fig 4.5 Maximum and minimum values of contact force between wheel and 

rail due to wheel flat (Neilsen and Igeland, 1995) 

(Contact forces are normalized with respect to half the static axle load. 

(a) one wheel flat hit midway between sleeper and (b) one wheel flat hit 
above sleeper. O, Wheel flat 60mm, A. wheelflat 90mm) 
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Fig 4.6 (a) Load history graph of 5 th span under 9 th and 10 th wheel 



Fig 4.6 (b) Load history graph of 10 th span under 9 th and 10 th wheel 
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Fig 4.7 Rail Response for Case 1 (a) - Rigid supports 

(a) node 10, (b) node 33 and (c) node 35 
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Displacement (m) 






Acceleration Response for Case 1 (a) - Rigid supports 
(a) node 10 , (b) node 33 and (c) node 35 






Fig 4.10 Acceleration Response for Case 1 (b) - Rigid supports 
(a) node 10, (b) node 33 and (c) node 35 







Displacement (m) 


x 10' 5 



45 





Displacement (m) 



0 0.5 1 1.5 2 2.5 3 3.54 



Fig 4.12 Rail Response for Case 2 (b) - Rigid Supports 

(a) node 1 0 , (b) node 33 and (c) node 35 
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(c) 

Fig 4.14 Rail Response for Case 2 (b) - Flexible supports 

(a) node 1 0, (b) node 33 and (c) node 35 
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Fig 4.15 Rail Response for Case 2 (b) - Rigid supports 
(a) node 82, (b) node 83 and (c) node 84 









(c) ZOOMED VIEW OF A FLAT WHEEL WITH NORMAL WHEEL 
Fig 4.16 Strain gauge Response obtained at test site 
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Fig 4. 1 7 Typical flat on a wheel 
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CHAPTER 5 


CONCLUSIONS 

A finite element model has been developed during the present work to describe the rail 
dynamics under a moving train. The model can incorporate train wheels with flats. The 
wheel flat has been modeled as a 'lift-hit' phenomenon. The rail is taken to be multi-span 
continuous beam and is modeled through 2-D beam elements. Natural frequencies and 
mode shapes are obtained and the results are validated by comparison with results 
obtained through analytical formulations. The train is modeled through multiple moving 
loads. The model is found to correctly represent the experimentally observed test 
patterns. Comparison of response patterns due to normal wheel loading and wheel-flat 
loading has been made. Distinct features have been found in the response patterns, which 
can aid in development of algorithms for wheel-flat detection. Displacement response is 
found to provide better information on the presence or absence of flatness in wheels. The 
acceleration plots appear to contain a large number of frequency components and need to 
be analyzed more rigorously, than done in the present work, for extraction of meaning 
information. 

The present work is restricted to a single flat on a single wheel on the train. The model 
developed during the present study can be readily employed to carry out further 
investigations on the rail response under the influence of multiple flats on the same 
wheel. Response under wheel-flats on more than one wheel also needs to be investigated. 
Such a study will help in designing better flat detection algorithms. Further work also 
needs to be carried out study the influence of various rail parameters like sleeper 
distance, stiffness and damping on the vibration magnitudes. 
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